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Abstract 

We propose a new model for explaining the observations of more than mass proportional heating 
of heavy ions in the polar solar corona. We point out that a large number of small scale intermit- 
tent shock waves can be present in the solar corona. The energization mechanism is, essentially, 
the ion reflection off supercritical quasi-perpendicular collisionless shocks in the corona and the 
subsequent acceleration by the motional electric field E = — (l/c)V x B. The acceleration due 
to E is perpendicular to the magnetic field, in agreement with observations, and is more than 
mass proportional with respect to protons, because the heavy ion orbit is mostly upstream of the 
quasi-perpendicular shock foot. The observed temperature ratios between O^^ ions and protons 
in the polar corona, and between a particles and protons in the solar wind are easily recovered. 
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The heating of the solar corona to temperatures of the order of 10^ K and more is one 
of the outstanding problems of solar physics. Beside the high temperatures, Soho/UVCS 
observations have shown that heavy ions in polar corona, like O^^ and Mg^+, are heated 
more than protons, and that heavy ion heating is more than mass proportional; further, 
the perpendicular temperatures Ti are much larger than parallel temperatures T\\ 

d of coronal holes. 
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J]. In addition. 



As a consequence of magnetic mirroring in the diverging magnetic fie 
heavy ions are observed to be faster than protons in the solar wind 
the coUisional coupling with protons up to 1.32 Rq indicates that the Mg^"*" heating has 
to be faster than minutes j^. These observations give stringent contraints on the coronal 
heating mechanism. Ion cyclotron heating has been considered since long (e.g., ^, 6, 7, 8|), 
but some details are not yet fully understood. The comprehension of coronal heating is of 
general physical interest, as the sun serves both as a huge plasma laboratory and as a model 
for a large class of stars. 

Shock waves are considered to be common in the chromosphere/transition region and 

nnnnn 

in the corona (e.g., [9|, llOj, lUl, ll2|, ll3[). For instance, photospheric convection leads to the 



emergence of small magnetic loops, which lead to magnetic reconnection with the network 
magnetic field; small scale plasma jets are formed in the reconnection regions, and fast 
shocks can form when jets encounter the ambient plasma 1^, [ll|, [l^. Indeed, recent X-ray 



Hinode and UV Stereo observations have shown that many more plasma jets are present 
in the polar corona than previously thought [15|, \\\ . Therefore, a large number of 
small scale, intermittent shocks can form in the reconnection regions and propagate toward 
the high altitude corona. Recent numerical simulations show that bursty, time de pen dent 
reconnection in solar flares can eject many plasmoids and create oblique shocks [l8[]. In 
the high corona, magnetic reconnection happens when current sheets form because of the 
evoM,, 00.0... .tn,ct.e. Q, wKUe .a,-e scale sKoC ~ tKe corona because 



solar flares and of the emergence of coronal mass ejections [20 



2l| . Such shocks are detected 



as type II radio bursts [22,]. In the case of solar flares, the associated reconnection outflow 



termination shocks can 
fraction of a second 14 



De so strong as to accelerate electrons to 100 keV energies in a 



23| . and in some cases both the upper and the lower termination 



shocks are identified in radio observations |24l ]. 

In the low /5, nearly coUisionless corona, a shock wave is formed when a superAlfvenic 
plasma flow having velocity Vi > Va collides with the ambient coronal plasma. Here, the 
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plasma (3 is given hy (3 = STip/B"^^ where p is the total plasma pressure, B is the magnetic 



field magnitude, Vi is the plasma velocity upstream of the shock, and Va = B j ^^-np is 
the Alfven velocity, with p the mass density. The Alfvenic Mach number is defined as 
Ma = Vi/Va- We notice that although the typical Alfven velocity in the corona, of the 
order of 1000 km/s, is larger than the observed jet velocities of 200-800 km/s 15|, |25| . 
the Alfven velocity in the reconnection region can be much lower, since B is weaker there. 
Indeed, the reconnection regions are characterized by current sheets, naagnetic field reversals, 
and magnetic quasi-neutral lines. For instance, Tsuneta and Naito IJ] argue that an oblique 
fast shock is naturally formed below the reconnection site in the corona, see their Figure 1. 
The plasma velocity in the reconnection outfiow region between the slow shocks is of the 
order of Va in the infiow region, that is much larger than Va in the outfiow region, and this 
leads to the formation of shocks (e.g., ay. 120|). _ 



Previously, shock heating of coronal heavy ions was considered by Lee and Wu 
but mostly in connection with subcritical shocks. Here, we propose that the more than 
mass proportional heating of heavy ions in polar coronal holes is due to ion refiection at 
supercritical quasi-perpendicular shocks and to the ion acceleration by the V x B electric field 
in the shock frame. In this connection, we notice that more than mass proportional heating 
of a particles and O^^ has been observed in the solar wind by the Ulysses spacecraft, between 



2.7 and 5.1 AU, downstream of interplanetary shocks, most of which were supercritical 
(see also Ref. |27|). 



26| 
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t is well known both from laboratory 28|, l29|] and from spacecraft experiments (e.g., 30, 



321]) that above a critical Mach number ^2.7 for perpendicular collisionless shocks 



(less than 2.7 if the shock is quasi-perpendicular), a fraction of ions, which grows with the 



Alfvenic Mach number 



33| , is refiected off the shock, leading to the so-called supercritical 



shocks. When the angle Obu between the shock normal (pointing in the upstream direction) 
and the upstream magnetic field is larger than about 45°, the refiected ions reenter the shock 
after gyrating in the upstream magnetic field. Such shocks are termed quasi-perpendicular. 
Conversely, for Obu < 45°, the refiected ions propagare upstream, forming the ion foreshock 
which characterizes the quasi-parallel shocks. The critical Mach number can decrease 
below 1.5 for oblique shocks in a warm plasma [34], so that ion refiection is a relatively 
common process. Ion refiection can be considered to be the main dissipation mechanism by 
which collisionless shocks convert the fiow directed energy into heat, while the electrons are 
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Upstream [ 



Downstream 




FIG. 1: Schematic of the magnetic field profile of a supercritical quasi-perpendicular collisionless 
shock. The main features like the magnetic foot, the ramp, and the magnetic overshoot are 
indicated. 



heated much less (typically, one tenth of proton heating) 31 1. 

For the solar corona, we consider a quasi-perpendicular supercritical collisionless shock, 
and we assume a simple one dimensional shock structure. The upstream quantities are 
indicated by the subscript 1, and the downstream quantities by the subscript 2. We adopt 
the Normal Incidence Frame (NIF) of reference, in which the shock is at rest, the upstream 
plasma velocity is along the x axis and perpendicular to the shock surface, Vi = {Vxi, 0, 0), 
the upstream magnetic field lays in the xz plane, Bi = (i?^,!, 0, -B^i), so that the motional 
electric field E = —V x B/c is in the y direction, Ey = VxiBzi/c. We further assume that 
Bzi ^ Bxi {Obu — 90°), in order to simplify the discussion. An order-of-magnitude estimate 
of the energy gained by ions after reflection at the shock can be obtained by approximating 
the reflected ion trajectory with a circle of radius r^, wit h rr, the ion Larmor radius, in 



the upstream magnetic field. Assuming specular reflection 30|, ImI, l35| , on average the ion 



velocity at the reflection point is perpendicular to the shock and along the x axis. Keeping 
in mind that ion reflection gives rise to a non adiabatic displacement in the y direction, the 
work W done by the electric field is 

W = q.EyAy, (1) 

where Ay ~ 2ri. For specularly reflected ions, the Larmor radius has to be evaluated with 
the upstream flow speed (neglecting the thermal velocity of the incoming ion distribution), 
i.e., v± ~ \Vxi\, so that 

VxiB^i rriiVxiC 



W = qiEy X 2rL = 2qi- 



q^B. 



(2) 



zl 



which yields W = 2miV^i. This estimate shows that the energy gain is mass proportionaL 
A more detailed calculation yields a more precise result, and shows that heavy ion heating 
is more than mass proportional. In order to do this, we remind that a distinctive feature 
of quasi-perpendicular collisionless shocks is the formation of a "foot" in the magnetic field 
profile in front of the main magnetic ramp, the latter culminating in the magnetic overshoot. 
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36| . see Figure 1. The 



beyond which the downstream values are gradually attained 
foot is due to the population of reflected and gyrating protons, which causes an increase in 
;he plasma density, and, because of the magnetized electrons, in the magnetic field strength 



291 . Even if the solar corona composition encompasses several ion species, the foot extent 
in the upstream direction is determined by the proton gyroradius, since protons are the 
major species. We define -Bfoot = (1 + b)Bzi, with b depending on the ion reflection rate; b 



can be estimated to be of the order 0.5-1 for typical shocks in the heliosphere [31|, |32|, |36 |. 
Direct observations in space show that very strong fluctuation levels are found in association 
with collisionless shocks. However, in what follows we will neglect fluctuations and we will 
consider only the average quantities, in order to set the stage. Taking into account the 
fact that the orbit in crossed electric and magnetic fields is a trochoid, we start from the 
equations of the particle trajectory. We assume the magnetic field to be along the z axis. 



and set the origin of coordinates at the point of ion reflection, with t = (e.g., Ref. 37l|): 



x(t) = -^sin((]t) + ^t (3) 

y(^) = ^[l-cos(^]^)] (4) 

where Q = qiB/rriiC, and B the local magnetic field. The corresponding particle velocity is 

cE 

v^{t) = -v±cos{nt) + —^ (5) 

B 

Vy(t) = v±sm{Qt). (6) 

Specular ion reflection implies that at t = the ion velocity is opposite to the incoming 
plasma velocity, Vx{t = 0) = —Vxi, whence 

v± = Vxi + — ^ = K.! + — = V^i——. (7) 

-Dfoot -Dfoot 1 + 

The reflected ions meet again the shock surface, at x = 0, at a later time ti > 0, correspond- 
ing to 

^sin(f]tO = ^t,. (8) 




FIG. 2: (Color online) Projection in the xy plane of the trajectories of hydrogen and oxygen ions 
reflected at the shock ramp. The motional electric field is also indicated. As in Figure 1, the 
vertical dashed lines separate the main magnetic field regions, such as the upstream region, the 
foot, the ramp, and the downstream region. 

Upon inserting the values of v± and of Ey in the above equation we obtain 

sin(nti) = (9) 

whose numerical inversion yields Vtti = 2.27885 for 6=1, ftti = 2.1253 for b = 0.5, and 
Qti = 1.8955 for 6 = (see below). At this time the particle will have moved in the y 
direction by an amount given by 

Ay(ti) = ^[1 - cos(^]^0] = ^^[1 - <^osint,)]-^±^^. (lo) 

il QiB^i (1 + by 

This displacement in the y direction determines the energy gained by reflected ions during 
the gyromotion in the field Ey-. 

W = q.EyAy = m,V^\[l - cos(fiti)]^^^ = ^^[1 " cos{nt,)]^m,V,\ (11) 

Taking into account the fact that protons move in the foot magnetic field -Bfoot, we can 
assume that b ^ 0.5-1. In such a case, 1 — cos(f2ti) = 1.65035 for 6=1 (or 1.52652 for 
b = 0.5), so that 

Wpc^^x 1.65035 X i^nipV^^) . (12) 
6 



On the other hand, for heavy ions hke O^^ most of the trajectory is upstream of the foot, 
see Fig. 2, in the unperturbed plasma where B B^i. Then we can set 6 = with good 
approximation, and obtain 1 — cos(i7ti) = 1.319, so that 

l^heavy ^ 4 X 1.319 X (^m.l^,^ ) . (13) 

Here we can see that, with respect to protons, heating is more than mass proportional For 
6=1, the ratio of the heavy ion energy gain over the proton energy gain is VFheavy/Wp — 
2.13 X mj/mp, while for h = 0.5 we find Wheavy/Wp ~ 1.55 x rrii/mp. Varying the value of 
b between 0.5 and 1 yields an 0^+ temperature about 25-34 times larger than the proton 
temperature (mo — IGrrip), in good agreement with Soho/UVCS observations which give 
To5+ /Tp = 27-37 [5I. Also, assuming a typical value of b = 0.5, we can easily recover the 
temperature ratios observed in the solar wind for helium {rria — 4:mp), where T^/ Tp — 6 is 



observed in those solar wind periods which are found to be the least collisional 



38j, thus 



reflecting more closely the coronal conditions. At the same time, for b = 1, T^/Tp ~ 8 is 



38|. 



obtained, a value which is also observed and reported in Figure 1 of Ref. 

On the other hand, heating is essentially perpendicular, since it is due to the motional 
electric field Ey which is perpendicular to the magnetic field by definition. This allows 
to understand the observed strong temperature anisotropy with Tj_ ^ Ty. In addition, a 
single shock encounter is required to accelerate the ions, and the acceleration time is on the 
scale of the ion gyroperiod, so that the heating mechanism is very fast, as required by the 
observations in Ref. [5|. 

Typically, the reflection rate for protons is found from numerical simulations to be 20- 
30%. The energy gained by reflected ions is distributed to the transmitted ions by wave 
particle interactions 31i , so that we can assume for the bulk of protons a heating rate about 
1/4 of Wp. Let us define the heating efficiency, rj, for protons, as the ratio of the energy gain 
of both reflected and trasmitted ions over the upstream thermal energy ^nipvff^i: 

,^J^^^f(f'-'fKo,2M^ (14) 

^rripV^i^i 4 X ^rripVf-^i 

where Mg = Vxi/vthi is the sonic Mach number. In the low (3 corona, the thermal speed is 
much less than the Alfven speed, so that Ms ^ Ma- However, in the reconnection outflow 
region the magnetic field is weaker than in the ambient corona, so that for a first estimate 
we assume that the sonic Mach number is of the same order of the Alfven Mach number. For 
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instance, assuming that Mg = 7, we can see that the efficiency for protons equals rj ^ 30. 
In other words, the crossing of two or three shocks might bring the chromospheric plasma 
from temperatures of the order of 10'^ K to coronal temperatures of the order of 10® K. As 
a general trend, we can say that the stronger the magnetic field in the reconnection inflow 
region, the larger the plasma velocity and the Mach numbers in the reconnection outflow 
region, and the larger the heating efficiency. On the other hand, several shock crossings may 
be required for the high altitude corona to reach the observed temperatures, and we notice 
that large scale shocks associated with coronal mass ejections and type II radio bursts can 
propagate all the way into polar corona. 

Clearly, the present model has to be further developed, since a wide range of different 
Mach numbers, plasma jS, shock normal angles dsn, fluctuation levels, and heating ratios 
can be envisaged in the corona. In addition, coUisional coupling with protons can decrease 
the obtained temperature ratios, as may be the case of Mg^"*". Further, multiple reflected 



43(1, a phenomenon which would increase the 



ions at a single shock can also be envisaged 
heating efficiency. 

Our model leads to the prediction that a fraction of heavy ions comparable to the proton 
fraction is also reflected at quasi-perpendicular shocks. In coUisionless shocks, an electro- 
stat ic p otential barrier A0 ~ mp{y^i — V^^)/2e arises, which slows down the incoming ions 



29|, 



39, 



44| . In simple discussions of ion reflection, ions are expected to undergo specular 



reflection if their kinetic energy in the shock frame is less than the potential energy barrier 
gjA0: however, for heavy ions this is usually found only for a tiny fraction of the upstream 
velocity distribution. Nevertheless, experimental evidence of a particle reflection at Earth's 
quasi-perpendicular bow shock has been reported by in Ref. [40|, while evidence of a par- 
ticle specular reflection off the quasi-parallel bow shock was reported in Ref. 41|. On the 
other hand, laboratory experiments show that increasing M^, the number of reflected ions 
increases while the potential jump decreases [29^, contrary to expectations if the potential 
jump would be the only cause of ion reflection. This shows that ions are not simply reflected 



by the average potential jump across the shock, and that also the magnetic 



fluctuating electric and magnetic overshoots play a role for ion reflection 



3, 



35, 



bot and the 



39 



42|. In- 



deed, cross shock electric fields measured by the Polar spacecraft at Earth's bow shock show 



that the potential A(/) is strongly spiky and fluctuating 



Further, quasi-perpendicular 



shocks also exhibit cyclic reformation, which implies time-depending electric and magnetic 
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overshoots and reflection rate 



33 



45|. Indeed, the Cluster spacecraft have recently shown 



that ion reflection is highly unsteady 45|], so that the strong variations in the shock structure 
can also induce heavy ion reflection. 

In conclusion, we have considered the heavy ion energization due to the ion reflection 
off quasi-perpendicular shocks. Fast, supercritical shocks are assumed to form because of 
reconnection of small scale magnetic loops at the base of coronal holes, like those associated 
with polar coronal jets, and because of the merging of magnetic structures in the higher 
corona. The energy stored in the coronal magnetic fleld is transformed to bulk kinetic 
energy by reconnection, and into heat and heavy ion heating by the quasi-perpendicular 
shocks which form in the reconnection outflow region. Our model can explain both coronal 
heating and the more than mass proportional heavy ion heating observed by Soho/UVCS. 
In addition, this heating mechanism is strictly perpendicular to the magnetic fleld and it is 
very fast (a single shock encounter is needed); most heating goes into the ions, with electrons 
undergoing an energy gain which is about an order of magnitude smaller than that of protons 
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32l |. Further, the strongly anisotropic heating with Tj_ ^ Ty can give rise to efficient 



ion cyclotron emission; this phenomenon is actually observed downstream of the Earth's 
bow shock. Indeed, recent Stereo, Helios, and Venus Express data show that ion cyclotron 
waves are probably present in the corona [46]. Therefore, quasi-perpendicular collisionless 
shocks can be the source of ion cyclotron waves in the corona, too. These waves later on can 
heat locally the solar wind by cyclotron resonance dissipation, as suggested by a number of 
observations. 
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